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Several recent exact diagonalization calculations have established that the Anderson-Hubbard 
model has a disorder-induced zero bias anomaly (ZBA) (also called a disorder-induced pseudogap) 
in the density of states. In order to understand the physics of the ZBA, we study a simplified 
problem — an ensemble of two-site molecules with random site energies — for which analytical results 
are possible. For this ensemble, we examine how the ZBA forms in both the weakly correlated (mean 
field) and strongly correlated limits. In the weakly correlated case, the ZBA can be understood as 
the result of level repulsion between bonding and antibonding molecular orbitals. A similar level 
repulsion occurs in the strongly correlated case too, but a larger contribution to the ZBA comes 
from the suppression of a triplet excitation mode. This inherently many-body mechanism does not 
have a counterpart in mean- field models. 



I. INTRODUCTION 

A number of recent papers have shown the exis- 
tence of a disorder-induced zero bias anomaly (ZBA) 
in the Anderson-Hubbard model (AHM) in one and 
two dimensions <i"— These calculations have revealed that 
there is a V-shaped dip in the density of states at the 
Fermi energy ep- This dip is produced by the response of 
the inelastic self-energy to the disorder potential. ^'^ Such 
a mechanism is well-understood in conventional metals 
and insulators, where the effect was explained at the level 
of Hartree-Fock theory by Altshuler and Aronov.^ How- 
ever, strong correlation effects are generally important in 
the AHM, and the Altshuler- Aronov mechanism is thus 
insufficient for this casci^ 

The AHM is the standard model for strongly- 
correlated systems with disorder. Like the Hubbard 
model, electrons are assumed to move on a tight-binding 
lattice of atomic-like orbitals. A zero-range intraorbital 
Coulomb interaction U is included, but longer range in- 
terorbital interactions are neglected. Strong correlations 
are important when the intersite hopping matrix element 
t is small relative to U . The AHM differs from the Hub- 
bard model by the addition of disorder, which is intro- 
duced by selecting the orbital energies from a random 
distribution of width A (the subscript i labels sites in the 
atomic lattice). The Hamiltonian is 

= -t ^ cJ^Cj^ + ^ {e,h, + U fii^fi,,,) (1) 

where {ij) restricts the sum to nearest neighbor sites, ni„ 
is the number operator for site i and spin a, hi = hi^, 
ande, e [-^A, iA]. (We use O to indicate the operator 
form of an observable O; thus rii — {hi).) In this model, 
the ensemble-averaged density is n = 1 (i.e. the band is 
half-filled) for £_f = f • 




FIG. 1: (color online) Comparison of densities of states 
for the Anderson-Hubbard model using different approxima- 
tions. Results are shown for the nonmagnetic Hartree-Fock 
(NMHF), unrestricted Hartree-Fock (UHF) and exact diago- 
nalization (ED) of small clusters. Hartree-Fock calculations 
are ensemble-averaged self-consistent calculations for 10 x 10 
lattices and 1000 impurity configurations. Exact diagonaliza- 
tion calculations are for a 12-site lattice and 1000 impurity 
configurations. Model parameters are A = 20, (7 = 8, t = 1 
and ef ~ U/2, corresponding to half-filling. 



The conventional Altshuler- Aronov theory predicts 
that the Hartree and exchange self-energies make posi- 
tive and negative contributions to the density of states 
a,t ep respectively/^ The exchange self-energy is typi- 
cally much larger than the Hartree self-energy, and the 
net result is a depletion of states aX ep- However, the 
AHM has a zero-range interaction for which the exchange 
self-energy vanishes. Altshuler- Aronov theory predicts a 
peak in this case, which is illustrated by the nonmag- 
netic Hartree-Fock calculations in Fig. [T] This is in 
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FIG. 2: (color online) Low energy excitation spectra as a func- 
tion of site energy for the degenerate two-site model. Lines 
represent peak position uj in the tunneling density of states, 
plotted as a function of ei for t = 0.5. (a) Excitation spec- 
trum for the mean-field Hamiltonian ((3| with V — 0.5 and 
£2 = ei; (b) Excitation spectrum for exact-diagonalisation 
calculations with 62 = ei — (7 and (c) £2 ~ ^i- Both (b) and 
(c) have U = 12. 



contrast to the V-shaped dip found in exact diagonal- 
ization calculations!^-"— The Altshuler-Aronov prediction 
assumes a nonmagnetic ground state, and a number of 
unrestricted Hartree-Fock calculations have found a V- 
shaped dip at ep^'^— in the magnetic phase<^ While the 
unrestricted Hartree-Fock results are qualitatively simi- 
lar to the exact diagonalization results, there are some 
important differences. Notably, the ZBA in the unre- 
stricted Hartree-Fock calculations grows with increasing 
U , eventually forming a broad soft gap when U is suffi- 
ciently large. In contrast, the ZBA in exact diagonaliza- 
tion calculations saturates for large U (provided U < A; 
a Mott gap opens for t/ > A), and empirically has a 
width ^ t.^'- Densities of states based on the different 
approximations are illustrated in Fig. [TJ 

We note that the above discussion ignores the low en- 
ergy soft gap^i^iiS that has been inferred from exact diag- 
onalization in one dimension, and found in unrestricted 
Hartree-Fock calculations in one and three dimensions. 
This gap appears on a scale |a; — EfI ^ 0(0. It), and has 
been ascribed to long range correlations.'^'^ The current 
work examines the two-site AHM where long range cor- 
relations are absent, and there is no soft gap. 

The advantage of the two-site AHM is that it is sim- 
ple enough that analytical results are possible, and yet 
is sufficiently rich to explain much of the physics of the 
ZBA in larger systems. Here, our goal is to compare 
the two-site AHM to a simple mean-field two-site model 
in order to answer the question: in what way is the ZBA 



in strongly correlated systems different from that in con- 
ventional metals? 

Our main results are summarized in Fig. [21 In 
Fig. [H^a) , we plot the low-energy excitation spectrum for 
a pair of sites with energies ei = £2 using the mean- 
field model described in Sec. [Hi This model is meant to 
illustrate the conventional Altshuler-Aronov mechanism 
for the ZBA. Without interactions, the hybridization of 
atomic orbitals leads to a level splitting of 2t between 
bonding and antibonding orbitals. With interactions, 
there is a range of site energies near sp where the level 
splitting is larger than 2t. This enhancement of level 
splitting (i.e. this level repulsion) shifts spectral weight 
away from ep and is the origin of the ZBA in this model. 
In Fig. [^fb), an alternative mechanism for shifting spec- 
tral weight away from e^? is presented. In this case, exact 
results for the low-energy excitation spectrum of the two- 
site AHM are shown. We have taken ei = £2 + U, which 
means that the lower Hubbard orbital of site 1 is degen- 
erate with the upper Hubbard orbital of site 2. Here, 
the spectrum has three excitation poles, the middle of 
which is a triplet excitation. The gap which is evident 
in the triplet spectrum is one of the main reasons for the 
pronounced ZBA in the two-site AHM, and is an inher- 
ently many-body mechanism that lies outside the mean- 
field Altshuler-Aronov paradigm. Finally, in Fig.[5Jc), we 
show that interactions have little effect on the spectrum 
if we consider the case of degenerate orbitals ei — ^2- In 
this case, excitation spectra are shifted by ~ 0{t^ /U) 
from their noninteracting values. 

The goal of this paper is to explore the physics behind 
these results. We discuss the mean-field mechanism for 
the ZBA in Sec.|TTl and emphasize in particular the role of 
level repulsion. We then derive, in Sec. lIIIi an expression 
for the ensemble-averaged density of states for the two- 
site AHM. Finally, we discuss in Sec. IIVI the different 
mechanisms by which the ZBA found in Sec. IIIII arises. 



II. ZBA IN MEAN-FIELD THEORY 

It is worth reviewing briefly how the ZBA arises in 
conventional metals. A variety of physical explanations 
for the Altshuler-Aronov ZBA have been given f^ii^iiS, and 
in this work we adopt the language of level repulsionJ^ 

We consider an ensemble of two-site AHMs with ran- 
domly chosen site energies. Since we restrict ourselves to 
nonmagnetic solutions of the Hartree-Fock equations, a 
V-shaped ZBA is possible only if a nonlocal interaction is 
included. We therefore add a repulsive interaction Vnih2 
to the Hamiltonian. In Hartree-Fock theory, 

Vhih2 V{nih2 + n2ni) 

~Vj2{{clc2a)clc^, + h.c). (2) 
a 

The first and second terms are the Hartree and exchange 
contributions respectively, and there is an additional 
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by the exchange self-energy. This enhanced level repul- 
sion, by itself, does not lead to a dip in the density of 
states; it is necessary that the amount of level repul- 
sion depend on the values of Ex,± relative to e^-. At 
zero temperature, Eq. ([5|) shows that t is different from 
t only if Ex,+ > sf > Ex-, as iUustrated schemat- 
ically in Fig. [3] This has the effect of pushing states 
away from ep , as shown numerically for the case ei = £2 
in Fig. IHa). In this language, the ZBA in conventional 
metals is understood as level repulsion between filled and 
empty molecular orbitals near ep- 



III. APPROXIMATE DIAGONALIZATION OF 
THE TWO-SITE AHM 



FIG. 3: (color online) Zero bias anomaly for self-consistent 
solutions of the mean-field Hamiltonian Q. (a) The density 
of states for t = 1 and diflferent A and V . (b) The exchange 
self-energy is nonzero for configurations of ei and £2 for which 
-B+ > ef > This leads to an enhanced level repulsion 

relative to configurations, as in (c) where E+ and E- are on 
the same side of ef and the exchange self-energy vanishes. 



Hartree contribution ^W^^nifii from the on-site inter- 
action. The Hartree contribution to the density of states 
is small for weak disorder^ but is central to the physics 
of the Coulomb gap for large disorder; the exchange con- 
tribution is largest for weak disorder, and underlies the 
Altshuler-Aronov mechanism for the ZBA. Because our 
goal is to contrast the Altshuler-Aronov mechanism with 
the physics of the the AHM, we discuss only the exchange 
term. 

Neglecting the Hartree contributions, we obtain the 
mean-field exchange Hamiltonian 



H 



X 



E 



t 



(3) 



where the renormalized hopping matrix element is t 
t -I- l^(c|g.C2o-). The eigenergies of Hx are 




Ex±= - , 

■ 2 V V 2 

and a straightforward calculation yields 

,t ^ ~!{Ex,^)-!{Ex,-) 



E 



X. 



Ex,- 



(4) 



(5) 



Equations @ and ([S]) allow i to be determined self- 
consistently for each (ei, 62) pair. The ensemble-averaged 
density of states for this model exhibits a ZBA, as shown 
in Fig. 13 

The term "level repulsion" refers to the fact that the 
level spacing between molecular eigenenergies is greater 
than the level spacing between the original atomic ener- 
gies, namely Ex.+ — Ex^- > |ei — £2]. For a repulsive 
interaction V, t > t and the level repulsion is enhanced 



We now turn to an approximate solution of the two- 
site AHM that preserves strong-correlation physics. We 
work in the strongly-correlated limit J7 3> t, where we can 
isolate terms that contribute to the density of states on 
the energy scale t. Higher order terms, which contribute 
on the scale t'^ /U, are neglected. We begin with a brief 
review of the atomic limit {t = 0), where interactions 
already have a nontrivial effect on the density of states, 
and then show how the density of states is modified by a 
nonzero t. 



A. Atomic Limit 

The density of states can be found exactly in the 
atomic limit t = 0. Each site is independent, and the 
ground state \Gi) for the ith site is 



|0), ep < £1 
= { It), ep-U <€.,< ep 
|2), e,<ep~U 



(6) 



We have assumed a weak Zeeman splitting so that spin- 
up states are preferred when there is an odd number of 
electrons. The spin-averaged retarded Green's function 
for the ith site is 



where cj^ 



2 ^ 

ma 

1 1 



|(m|c,,|G,)|2 



Eg. 
1^ 



E„ 



u+ -f Eg, - E„ 

(7) 

|m) are a complete 



CJ -I- iO, Ui = Y.ai'^i' 

set of excited states with energies Em, and E^ is the 
ground state energy. The spin-averaged density of states 
at site i is thus 



1 



Im GM 



(8) 



I-^)<5(c.-£,) + ^<5(a.-£,-C/). (9) 



This equation shows that (i) strong correlations split the 
local spectrum at each site into a pair of poles at and 
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ei + U and (ii) the weight of each pole depends on the elec- 
tron density at that site. We refer to the poles at and 
ei + U as the lower Hubbard orbitals (LHO) and upper 
Hubbard orbitals (UHO) respectively. It is worth empha- 
sizing that the energies of the LHO and UHO determine 
the total charge density at each site. From Eq. ([6]), 

0, Ef < (LHO and UHO above ep) 

1, ei < ef < £i + U (LHO below; UHO above) 

2, e; -f t/ < £F (LHO and UHO below ef) 

(10) 

At half- filling {ef = U/2), the ensemble-averaged den- 
sity of states is 
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-U -UJ 



(11) 



where Q{x) is the step function. The result (ITl]) is illus- 
trated in Fig. m This figure explicitly shows the spectral 
weight contributed by the LHO and UHO in their dif- 
ferent filling states. For this work, the most important 
aspects of the figure are (i) that both LHO and UHO con- 
tribute spectral weight at Ef for Ef ^ [U — 5 A, 5A], and 
(ii) that for this range of there is a "central plateau" 
where interactions enhance p{ef) relative to the nonin- 
teracting value A~^. 



B. Two-Site Case 

1. Preliminary Discussion 

The results of exact numerical calculations of the den- 
sity of states are shown in Fig. [5] for two cases: C/ < A 
and [/ > A. We track the evolution of the density of 
states as a function of in both cases. When U < A, 
there is a broad ZBA centered at ef for ef = 3 and 
Ep = 5. However, the ZBA is unresolvable when Ef is 
outside the central plateau. When U > A, there is a 
Mott gap at half-filling, and a ZBA forms as one dopes 
away from half-filling. This ZBA is qualitatively differ- 
ent from that found near half-filling for J7 < A and has 
a width of order /U . In this section, we focus on the 
large ZBA that appears near half-filling for U < A. 

The approach we take is to calculate the density of 
states p^^^i^iu)) for a single two-site AHM with site en- 
ergies ei and £2. The density of states is then averaged 
over all possible configurations. 



p{uj) 



iA 



dei 



-iA 



de2 P£i,£2 {^)- 



(12) 



-iA 



To simplify the analytic calculations, it is useful to 
partition the integration range [— ^A, ^A] into subranges 
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FIG. 4: Density of states in the atomic limit. The figure shows 
different contributions to the ensemble- averaged density of 
states in the limit t = for f/ = 12 and A = 20. The spectral 
weights contributed by lower Hubbard orbitals (LHO) and 
upper Hubbard orbitals (UHO) in their different filling states 
are shown. For comparison, the noninteracting density of 
states is A"'^ for — |A < lu < |A. The density of states 
in the "central plateau" is §A^^ and is thus enhanced by 
interactions, relative to the noninteracting case. The central 
plateau extends over \U — ^A, ^A] and is the region where 
the LHO and UHO coexist. 



U=6, A=12 



U=16, A=12 




FIG. 5: (color online) Exact numerical solution for the density 
of states. Density of states for (a) (7 < A and (b) [/ > A for 
different values of ef. In both panels, half-filling corresponds 
to Ef ~ U/2. Curves are offset for clarity. Densities of states 
are averaged over 10^ random disorder configurations. 



A= [-^A,0] and S = [0, iA], as illustrated in Fig. H] 
Sites whose UHO lies near ef belong to region A, while 
sites whose LHO lies near ef belong to region B. We 
have argued^^ that the ZBA comes from level repulsion 
between LHO and UHO on neighboring sites, and it is 
indeed suggested by Fig.^b) and (c) that the important 
configurations have ci £ A, £2 & B or ei £ B, 62 & A. 
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The simplest approximation is to treat these configura- 
tions carefully, while treating the other configurations in 
the atomic limit. As we show, this turns out to be suffi- 
cient to understand the essential physics of the ZBA. 

We denote by pxyi^^) the density of states ensemble- 
averaged over sites with ti £ X and €2 G 3^, 



(a) 



pxyiuj)=lm dei / de2 Pei.ea ('^)- (13) 
Jx Jy 

Paa{'-^) a-nd /9j3e('^) sue evaluated in the atomic limit, 
using Eq. (0), 



_ 1 

2A 



+ -eiU-\u;-eF\) 



(14) 



For w and ep near U/2 (half-filling), PAA+BBii-^) — 
|A~^. Using Pba{^) — Pab{^)^ the total density of 
states is 



(15) 



A more careful derivation oi p^j(^gis{uj) finds corrections 
to the atomic limit approximation on the energy scale 
\uj-eF\ < 0{t^/U). 

The next step is to evaluate 



Pba(.^^) 



— \pr [ dei [ (ie2lniG'ei,e2(a;), 

Jb J a 



with Gei,e2(w) the retarded Green's function averaged 
over sites and spins. It will be convenient to change in- 
tegration variables to 



e2 + U + ei 



e2 + U - ei 



Sf 



(16a) 
(16b) 



and write 
Pba = 



dy 



-A 



dxlm Gx,y{i^) (17) 



where the factor of 2 is the Jacobian for the transforma- 
tion, and the integration limits are 



A + A 



A- A 



y 



and 



A = 



A-U 



-I- 



The Fermi energy is written 



eF = ^ + 6. 



(18) 



(19) 



(20) 
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FIG. 6: (color online) Phase diagram of the two-site AHM. (a) 
Filling states as a function of ei and 62 for A = 20, U = 12, 
t — 1, and ef ~ U/2, corresponding to half-filling on average. 
Also shown are the regions A and B. Solid black lines indicate 
phase boundaries for t nonzero; dotted red lines at ef and 
Ef — U indicate phase boundaries for t — 0; dashed blue 
lines indicate boundaries of regions A and B. (b) Integration 



region BA for ef 
coordinates. 



ill + S, shown in terms of transformed 



This equation defines 6 such that half-filling corresponds 
to S = 0. Figure [6] illustrates the integration region and 
gives the graphical meaning of A, A, and S. 

The phase diagram Fig. |6] shows that there are three 
filling states in BA, with iV = 1,2, or 3 electrons. We 
now find the ground state wavcf unctions, energies, and 
phase boundaries for the different filling states. 



2. Ground states in region BA 

We will determine the ground state wavefunction in 
the region BA using a truncated basis set that discards 
high energy states. These high energy states modify the 
ground state wavefunctions and energies by 0{t'^ /U), and 
our approximation is consequently valid for U ^ t. 

In BA, the one-electron ground state in the atomic 
limit is |0 t) because ci > €2- Making t nonzero mixes 
in a small amount of | t 0), proportional to i^/(ei — £2)- 
However, in BA, ei ~ e2 -|- U, so the mixing is of order 
t'^/U and is neglected in our approximation. The one- 
electron ground state is thus 



|Gl)«|0t). 



(21) 



Similarly, the three-electron ground state is \G3) « | t 2). 

The two electron ground state is found by diagonaliz- 
ing the AHM Hamiltonian in the reduced basis {\s), |02)}, 
where 



1 
72 



(I n) - 1 ;t)) 



(22) 



is the singlet state. For ci ^ €2 + U ^ Ef, each of these 
basis states has an energy ^ 2£f — U . The discarded 
basis state |20) has an energy ^ 2£f + U , and the amount 
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TABLE I: Approximate A'^-electron ground states \GN) and 
their energies Eg for the two site model with A > U ^ t in 
the region BA. Variables x and y are defined in (|16p , and E- 
is defined in (l24l). 



N 


Eg 


-epN 


|G7V> 


Ground state when 


1 




- Ef 


lot) 




^ > \Jy'^ + 


2 




- Ief 


ai_ s) + 02- 


-102) 


|x| < Vy' + 2i2 


3 


ei + 262 




It 2) 







where x and ?/ are defined in and w = w — Si? . 
The ground state has three electrons for 



Xy - E <x < -y/y'^ + 2*2, 



(30) 



where the upper limit is the phase boundary between 
two- and three-electron states (c.f. TablelH, and the lower 
limit [c.f. Eq. ([T5|) ] is the boundary of region BA. Then 



of 1 20) mixed into the ground state by t is therefore ~ 
0{t'^/U), which we ignore. The Hamiltonian matrix in 
the reduced basis is 



H 



€l + €2 -V2t 
-V2t 2e2 + U 



(23) 



which has eigenenergies 



± 



£2 + C/ - £1 



and eigenstates 



|±) =ai±|s) +a2±|02) 



a 



i± 



^2± 



1 - a 



+ 2*2 (24) 

(25) 
(26) 



i± 



where y is defined in p^ . The two-electron ground state 
is \G2) = |— ). The different ground states, their energies, 
and the phase boundaries between them are tabulated in 
TablelH The next step is to calculate the density of states 
for each filling state. 



3. Density of states for the three- electron ground state 

First, we calculate the contribution to the density of 
states from the 3-electron ground state. Throughout this 
work, we keep only terms with poles near ep, meaning 
that terms with poles near ei + U or £2 are discarded. 
The spin- and site-averaged Green's function is then 



|(±|cn|t2)|2 + |(±|c2tlt2)|^ 
^ ^ a;+ - (£1 + 2£2 + U) + E± 

|(*|c2tlt2)P + |(tt|c2ilt2)P 

UJ+ - (£1 + 2£2 + U) + (£1 + £2) 



,(27) 



where 



(28) 



\t) = ^ (in) + 1 ;t)) 

is a triplet state. Using a2± = 1— we reduce Eq. ([27] 
to 



1 



^ Q+ -x± ^2/2 ^ 2*2 8(2}+ -x~y 

(29) 



Pba 



2A2 



dy 



dx 



(1 - al+/2)5{Co -x+ yfy-^ + 2*2) 



+(1 - a2_/2)(5(w ~x- 7^2 + 2*2) 



+ -5{[b -x-y) 



(31) 



The integration over x is straightforward because of the 
delta-functions, which introduce the constraints uj < Q 
and 



, , l^l / 8*2 

\ + 6-\Cj\ K + &-\Cj\ 

2 2 ' 

\Cj\ *2 \ k + 6-\Co\ 

T-M 2 ' 



first term, 
second term, 
third term. 



The result for the first term is valid forO>£l'>— A — 5, 
i.e. for uj < ef and in the central plateau. In deriving 
these results, we have neglected terms of order *2/A and 
*2/A. We now integrate over y using 



dy 1-^ =7 32/±Vy^T2^ 



(32) 



to get the three electron contribution to the density of 
states, 



8A2 

+ e 



2A + A , , , / / 8*2 ' 

_ 2*' 

~ K + 5 



A . 2*2 
A-f <5- — 
\lo\ 



(33) 



To simplify the final expression, we have taken 
VA2 + 2*2 ss A and VA2 + 2*2 w A. 



4- Density of states for the one-electron ground state 

The derivation of the one-electron contribution to the 
density of states parallels that of the three-electron con- 
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tribution. The Green's function is 

1 r^i(±i4,|ot)p + i(±i4,iot)p 



UJ+ +62- E± 



\{t\ci\ot)? + m\ci,\ot)? 



UJ+ + €2 — (ci + £2) 

1 - af±/2 



+ - x T vV + 2t2 8 w+ - x + y 

(34) 



given by and 



PBA - OA 2 



. „ , rfv / dx 



(1 - a?+/2),5(cD - a; - Vj/' + ^f^) 



+ (1 - a?_/2)J((i - X + + 2t2) + - a; + ?/) 

Letting x — ?> — x, this is the same as p'^^i^^) for S — ^ —(5. 
Thus 



2^2 



A r ; > — 

A- (5/ V w 



(35) 



5. Density of states for the two-electron ground state 
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The integration region is Xy — S > x > ^/y"^ + 2t'^ with 
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FIG. 7: (color online) Density of states for the two-site 
ensemble-averaged Anderson-Hubbard model for different val- 
ues of t and 5. Results are shown for exact numerical solu- 
tion of the AHM (dashed black curves) and for the approx- 
imate result, Eq. (|40|) (solid red curves). Model parameters 
are 17 = 12, A = 20. Recall that ef = U/2 + 5 and that 5 = 
corresponds to half-filling. 



gives 



4A2 



2A + A 



8^2 

1-Re \/l - — ^ 



(39) 



Finally, the Green's function for the two-electron 
ground state is 



G 



2e ^ 1 r|(0i|cit|G2)|2-f |(0i|c2t|G2)|2 



Then, 



2 [ uj+^E^+e2 
|(t2|clJG2)|2 + |(t2|4JG2)|2 



± 



0J+ +E^- (ei + 262 + U) 
l-at/2 



2 ^ w+ - a; ± A/y2 + 2l2 



1 /-A-l^l r^2 / Q,2 



(36) 
(37) 



(38) 



where Xi = — min(a;j, -I- 6, y^y-^ + 2t'^), X2 = min(a;y — 
6, + 2t2). Perfo rming the integrations over x and y 



6. Total density of states 

Putting the results of the different calculations to- 
gether, we arrive at our final result for the total density 
of states (valid in the central plateau) 

p{oj) = P^^+bb('^) +2p^°^(w) -H24<=^(w) -H24°^(w) 
_ 3 3 \ U 6s 



+- 



1 - Re a/1 - 



8^2 



+ e |c:;| 



4*2 



U-2S, 



U -26, 



(40) 



where 6s = 6 sgn(a)) and Cj — uj ~ sf- The parameters A, 
A and 6 are defined in Eqs. (HH) and (pO| . 

As a check of Eq. (|^n)l . we let t 0, in which case 



p{uj) 



2A 
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which agrees with the previous atomic hmit calculation 
for the central plateau. Equation is plotted in Fig. [7] 
in comparison with the exact density of states determined 
from numerically diagonalizing the AHM. 

IV. DISCUSSION 

Figure [7] compares Eq. pOI) to exact disorder- averaged 
numerical calculations for the density of states. The the- 
ory works well for \uj\ > At^/U, out to the edges of the 
central plateau where it breaks down (for example, near 
uj = ±4 in Fig. EJa)). The theory neglects terms of or- 
der t'^ /U, and therefore fits the numerics better when t 
is smaller, as shown in Fig. [TJc) and (d). The fit for 
\u}\ < it'^/U is not especially good, but can be improved 
significantly by considering corrections of order t'^ /U that 
were neglected in the previous section; we haven't in- 
cluded these corrections because they complicate p{E) 
significantly without adding physical insight. The focus 
of this discussion is therefore |w| > it'^/U. 

The main qualitative idea that we emphasize in this 
section is that there are two distinct physical mechanisms 
that lead to the ZBA in Eq. (PH)) . Both mechanisms 
occur for configurations where the LHO of one site and 
the UHO of the other site are nearly degenerate with ep, 
namely for ei ~ €2 -I- J7 ~ £f or £2 ~ ei -I- C/ ^ ep- 
The first mechanism is similar to that outlined in the 
mean-field calculation in Sec. |TT| level repulsion, caused 
by hybridization of many-body states, shifts states away 
from £f- The second mechanism does not have a mean- 
field counterpart: level repulsion gaps the spectrum of 
low energy triplet excitations. 

The first mechanism underlies the second last term in 
Eq. go]), 

3|w| / I ^\ , , 

This term rises linearly from a) = and is peaked at 
cj = ±2\/2t, which defines the width of the ZBA in 
Fig. [71 In our calculations, this term comes from transi- 
tions between two-electron singlet states and states with 
one or three electrons, and it is the level repulsion be- 
tween the two-electron states that causes the ZBA. In the 
case, for example, where t\ and e2 -I- J7 lie near ep, there 
are two nearly-degenerate two-electron singlets, |s) and 
1 02); these hybridize as a result of the matrix element t 
to form bonding and antibonding many-body states with 
energies (from Eq. 

E±^2eF-V ± V2t. (42) 

Thus, the level repulsion between \s) and |02) shifts the 
many-body orbital energies up or down by 0{t). Starting 
from the two-electron ground state, with energy , one 
has transitions 

ai_|s)-t-a2-|02)''"4^' |a2), .43) 
ai-|s) +^2-102) "'"4'' \0a). 



We showed in Sec. IIIIBI that the three-electron energy is 
ei -|- 2e2 + U (which is approximately Sep — U), and the 
one-electron energy is £2 (approximately sp — U), so that 
the transition energies in Eq. are 

w±w£F±%/2t. (44) 

Because uj± are shifted by 0{t) away from sp, the den- 
sity of states at ep is reduced as t increases. As indi- 
cated above, this mechanism for depleting the low energy 
density of states is similar to the mean-field mechanism 
discussed in Sec.|lTl where level repulsion between molec- 
ular states on opposite sides of Ef increases the energy 
required to add or remove an electron. In this sense, the 
second-last term in (PU)) is Altshuler-Aronov-like. 

The second mechanism does not have a mean-field 
counterpart, and results in the last term in Eq. ()40|) 

This term varies as |£D|~"'^ down to the low energy cutoff 
at ~ At'^/U, and makes the dominant contribution to 
the shape of the ZBA. The cutoff comes from the bound- 
ary between the region BA and the region BB in Fig. IH 
where the approximate one- and two-electron wavefunc- 
tions used in deriving p{E) cease to be valid. 

In our calculations, Eq. ()45|) comes from transitions 
between one- or three-electron ground states, and two- 
electron triplet excitations. For the three-electron ground 
state, for example, these excitations have the form 

|t2)'^^|n), lt2)'A^|i). (46) 

As mentioned above, the three-electron energies are 
nearly independent of t; the triplet energies are also in- 
dependent of t, however, so that the transition energies 
are not shifted by level repulsion. The mechanism for 
depleting the low energy density of states in this term is 
therefore not that of Altshuler and Aronov. 

Instead, it is the fact that a gap in the triplet spec- 
trum opens as t increases that causes a depletion of 
states near ep (this gap was illustrated in Fig. [2]). This 
gap occurs for configurations of (ei,e2) that have one- 
or three-electron ground states in the atomic limit, but 
two-electron ground states when t is nonzero. For exam- 
ple, when €1 and €2 + U both lie slightly below ep, the 
atomic-limit ground state has three electrons and triplet 
excitations as in (|46| are possible. When t is nonzero, the 
two-electron ground state energy is reduced by 0{t), 
while the three electron ground state energy is reduced by 
0{t^ /U). For sufficiently large i, the two-electron ground 
state has the lower energy and the triplet excitation is 
eliminated (i.e. the only possible tunneling processes are 
to one- or three-electron final states). In summary, the 
ZBA in the final term of Eq. ([Ml occurs because the 
phase space for low energy triplet excitations is reduced 
when t increases. 

As we discus elsewhere^iii this calculation sheds light 
on the empirical observation, made for larger systems. 
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that the width of the ZBA is of order ti^ A nai've 
argument based on the disorder-free Hubbard model 
would suggest that the ZBA might have a conventional 
Altshuler-Aronov form, but with an effective exchange 
interaction Foff — 4:t'^/U, so that the ZBA should grow 
with increasing P/U. As we have said above, there are 
contributions to the density of states of this type; how- 
ever, we have just shown that a much larger effect, of 
order t, comes from configurations with the LHO and 
UHO on neighboring sites degenerate. 

We note that this explanation appears to contradict 
numerical evidence from the work of Chiesa et ali^ on 
two-dimensional clusters that a large ZBA persists far 
from half-filling and for large U, since configurations with 
a degenerate LHO and UHO do not occur in these cases; 
as we have shown in Fig. [Sj the ZBA vanishes rapidly 
(with increasing disorder) in the two-site model when 
Ei? is outside the central plateau. To check this, we have 
performed preliminary exact diagonalization calculations 
for larger clusters (up to 12 sites). These calculations find 
that the width of the ZBA is not linear in t when ep is 
outside the central plateau, and suggest that the physics 
of the ZBA changes far from half-filling. A more detailed 
study of how the ZBA evolves with doping needs to be 
undertaken. 



V. CONCLUSIONS 

In summary, we have found that the zero bias anomaly 
in the two-site Anderson-Hubbard model is the result of 



strong orbital hybridization in the two-electron ground 
state for configurations with ei ^ £2 + U ^ ep or with 
€2 ^ £1 + U ^ ep- Unlike in the conventional Hubbard 
model, this hybridization is not suppressed by the on- 
site interaction U, and leads to a level repulsion between 
molecular orbital energies of order t, rather than /U . 

The mechanism for the suppression of the tunneling 
density of states is, at least in part, different from in con- 
ventional mean-field models of interacting electrons. In 
mean-field theories, interactions cause a shift of molecu- 
lar orbital energies away from ep that leads directly to an 
increase in the energy required to remove or add an elec- 
tron. This also occurs in the Anderson-Hubbard model; 
however, there is an additional depletion of low energy 
spectral weight because the low energy triplet excitation 
spectrum is gapped as a result of orbital hybridization. 
This mechanism is physically different from that of Alt- 
shuler and Aronov. 
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